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A systematic analysis of the Burgers — Kardar-Parisi-Zhang equation in d + 1 dimensions by dynamic renormalization group 
theory is described. The fixed points and exponents are calculated to two-loop order. We use the dimensional regularization 
scheme, carefully keeping the full d dependence originating from the angular parts of the loop integrals. For dimensions less 
than d c — 2 we find a strong-coupling fixed point, which diverges at d = 2, indicating that there is non-perturbative strong- 
coupling behavior for all d > 2. At d = 1 our method yields the identical fixed point as in the one-loop approximation, and 
the two-loop contributions to the scaling functions are non-singular. For d > 2 dimensions, there is no finite strong-coupling 
fixed point. In the framework of a 2 + e expansion, we find the dynamic exponent corresponding to the unstable fixed point, 
which describes the non-equilibrium roughening transition, to be z = 2 + 0{e' i ), in agreement with a recent scaling argument 
by Doty and Kosterlitz. Similarly, our result for the correlation length exponent at the transition is 1/v — e + 0(e 3 ). For the 
smooth phase, some aspects of the crossover from Gaussian to critical behavior are discussed. 

PACS numbers: 05.40.+j,64.60.Ht,05.70.Ln,68.35.Fx 



I. INTRODUCTION 

In recent years there has been much theoretical interest 
in the dynamics of growing interfaces in random media 
pip) . This problem has received particular attention not 
only because of its technological importance and appli- 
cations but also as the simplest non-trivial example 
of dynamic scale invariance in a nonequilibrium system. 
Kardar, Parisi, and Zhang (KPZ) have proposed a sim- 
ple nonlinear Langevin equation |?| which has become a 
widely accepted and by now well-known description of 
the macroscopic aspects of a certain class of growth pro- 
cesses. This equation is closely related to a large vari- 
ety of other problems ranging from randomly-stirred flu- 
ids [Q (Burgers equation), dissipative transport J5| (the 
drivcn-diffusion equation), flame-front propagation ft 
(Kuramoto-Sivashinski equation), polymer physics 
the dynamics of a sine-Gordon chain ||, and the be- 
havior of magnetic flux lines in superconductors |lCj |. 

Through a simple transformation the Burgers-KPZ 
equation can be mapped onto the statistical mechanics 
of directed polymers in a random medium |8[p"H]. This 
problem seems to embody many of the features of more 
complex random systems such as spin glasses ]T^ ]. 

Due to the simplicity of its form, accompanied by a 
stunning complexity of its behavior, it has gradually 
taken on the role of an "Ising model of nonequilibrium dy- 
namics" . Therefore, any advances in understanding the 
behavior of the Burgers-KPZ equation may have broad 



impacts in both the fields of nonequilibrium dynamics as 
well as disordered systems. 

Assuming that a coarse-grained interface of a d— 
dimensional substrate may be described by a height func- 
tion h(x,t) with x € lZ d , KPZ proposed that the follow- 
ing nonlinear Langevin equation S , 



(1.1) 



governs the macroscopic (large-distance, long-time) dy- 
namics of a s toch astically growing interface. The first 
term in Eq. (1.1) represents the surface tension which 
prefers a smooth surface. The second term describes the 
tendency of the surface to locally grow normal to itself, 
and is entirely nonequilibrium in origin. The last term is 
a Langevin noise to mi mic the stochastic nature of any 
growth process. In Eq. (1.1), the average growth velocity 
has been subtracted so that the noise has zero mean, i.e., 
(?/(x, t)) — 0. In the simplest case, the stochasticity is 
then described by an uncorrelated Gaussian noise with 
the second moment 



<7j(x, t)r)(x', t')) = 2D6 d (x - x')S{t - t 1 ) . 



(1.2) 



where D characterizes the noise amplitude. 

The shape of the steady-state surface profile is com- 
pletely characterized by the steady-state distribution 
function 'Pf/i^x, £)]. The leading (second) moment of V 
is the truncated two-point correlation function 



C(x, t) = ([/i(x + x, to + *) - h(xo, to)] 2 



(1.3) 
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In Eq. ( |l.3| ) i the mean value of h has been implicitly sub- 
tracted, i.e., we use a moving reference frame suc h th at 
(d t h(x,t)) — 0. Because the equation of motion (1.1) is 
scale invariant ||, the correlation function assumes the 
scaling form 



C(x,t) = x 2x F(t/x z ), 



(1.4) 



where x = |x|. In Eq. (1.4), the roughness exponent x 
describes the scaling of the width of the interface, and the 
dynamic exponent z characterizes the spread in time of 
disturbances on the surface. In the limits y = and y — > 
oo, the scaling function F(y) becomes A — const, and 
By 2x / Z , (with B = const.), respectively, thus yielding the 
usual asymptotic scaling form of the correlation function, 



C(x,t = 0) = Ax 2x , 
C(x = 0, t) = Bt 2 *' z . 



(1.5) 



In cases where the nonlinear term in Eq. (1.1) is not 
relevant, one recovers the linear equation of Edwards and 
Wilkinson ]T^| . There the exponents are known exactly, 
Xq = (2 — d)/2, and zo = 2, as expected from the ensuing 
simple diffusion equation. 

In the more interesting case where the nonlinear term 
is important, one may derive the exponent identity 



X + z 



(1.6) 



which follows from the invariance of (|1T|) with respect to 
an infinitesimal tilting of the surface ||3| (h — ► h + v • x, 
x — * x — Avt). This invariance is known as "Galilean 
invariance" due to the corresponding invariance of the 
Burgers equation 0|. Therefore, there is at most one 
independent exponent to be determined. However, that 
last link needed to determine the exponents has proven 
to be quite elusive despite considerable effort. 

In the special case of d = 1, which is realized for exam- 
ple in the expansion of a domain boundary separatin g the 
two phases of a two-dimensional Ising model, Eq. ( |l.l| ) 
also satisfies a fluctuation-dissipation theorem (FDT). 
Then one can show |l]] that the stationary distribution 
function is 



P st [h(x)] = exp ■ 



v 

2D 



dx 



dh 

dx 



(1.7) 



Hence x — 1/2, as if the nonlinearity were absent. Given 
the relevance of nonlinearity in d = 1, one immediately 
obtains z = 3/2 through the exponent identity, Eq. (1.6). 

For dimensions d > 2 an interface described by 
the Burgers-KPZ equation is expected, on the basis 
of renormalization-group arguments to undergo a 

transition from an asymptotically smooth to a rough pro- 
file as the effective strength of the nonlinear coupling 
A 2 D/v 3 is increased. The exponents in the rough phase 
are only known from numerical simulations. The most re- 
cent values Hi seem to settle between the Wolf-Kertesz 



conjecture [15] x/ z = l/(2d + 1) and the one by Kim and 
Kosterlitz [16] x/ z = l/(d + 2). Despite those results 
from numerical simulations, there has been very little 
progress towards an analytical theory which goes beyond 
the original paper ra, at least on the basis of (refined) 
perturbational approaches. The main obstacle within a 
systematic treatment originates from a renormalization 
group (RG) flow towards a strong-coupling fixed point. 
Therefore, dynamic renormalization group (DRG) meth- 
ods have had very limited success. The main qualitative 
results obtained ||]4|] are (i) the interface is always rough 
(i.e., the nonlinearity is relevant) for d < 2 p7j , and (ii) 
there is a dynamic phase transition between a smooth 
phase (nonlinearity irrelevant, x — Xo,z = zq) and a 
rough phase where the nonlinearity is again relevant for 
d > 2. Quantitatively, DRG provides a perturbative ac- 
cess to the dynamic phase transition point via a d = 2 + e 
expansion. Unfortunately, DRG tells us virtually noth- 
ing about the most interesting aspect, namely the rough 
phase of the interface. The method breaks down be- 
cause the DRG fixed point describing the rough phase is 
a strong-coupling fixed point generally not accessible by 
any perturbational means. This problem persists even in 
d = 1 where the exponents themselves are already known 
exactly. 

Recently, there has been some progress in understand- 
ing the behavior at the roughening transition itself, which 
is described by an unstable fixed point in the DRG the- 
ory. Doty and Kosterlitz jll| have argued that the dy- 
namic exponent z c at the KPZ roughening transition 
equals 2 for all dimensions d > 2. Their argument is 
based on the mapping of the Burgers-KPZ equation to 
the equilibrium model of directed polymers in a random 
medium, and a standard scaling argument usually ap- 
plied for glassy system, also very similar to the deriva- 
tion of hyperscaling at critical points. Basically, they 
conclude that at a finite-temperature fixed point, one 
expects the scale of the free energy (which corresponds 
to the height in the Burgers-KPZ equation) to be given 
by the temperature, which is finite. This implies for the 
roughness exponent Xc = and z c — 2 via the exponent 
identity Eq. (1.6). This view is indeed supported by nu- 
merical simulations by Tang, Nattermann, and Forrest 
p9| . Based on the one-loop approximation, Nattermann 
and Tang Q] also discuss the interesting crossover scal- 
ing behavior of the Burgers-KPZ equation. 

The purpose of this paper is to go beyond the one- 
loop perturbation theory which is plagued by some arte- 
facts, and perform a two-loop calculation, in order to 
possibly clarify which aspects are generally accessible to 
a DRG treatment, and which are not, and also to es- 
tablish certain trends, which might provide clues even 
beyond the two-loop approximation. We shall apply the 
field-theoretical version of the DRG in the formulation of 
Bausch, Janssen, and Wagner 121122]. Using the dimen- 
sional regularization procedure 122 [for a massless theory 
such as provided by Eq. (1.1), requires a clear distinction 
between infrared and ultraviolet singularities, only the 
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latter to be included in the renormalization constants. 
Also, in order to be able to correctly describe the ex- 
actly known behavior in d = 1 dimensions, performing 
an e expansion turns out to be a very delicate procedure. 
For d < 2 our results indicate that such an expansion 
becomes inconsistent, the strong-coupling fixed point di- 
verging for d — > 2, which is very drastically different to 
the usual situations, as, e.g., in the paradigmatic <j> the- 
ory ^4j,^5|. On the other hand, for d > 2 the now un- 
stable fixed point describing the roughening transition 
may be treated in the framework of a 2 + e expansion. 
The situation here can in fact be compared to the similar 
treatment of the Oin) nonlinear a model ^6|, where the 
nonlinear coupling corresponds to the transition temper- 
ature, and is of order e, too. 

During the process of completing this paper, we be- 
came aware of a similar field-theoretic two-loop analysis 
of the Burgers-KPZ equation by Sun and Plischke p7| . 
These authors, however, apply a somewhat different pro- 
cedure, and some of their results are clearly incompatible 
with ours, as shall be discussed below. 

The outline of this paper is as follows. In Section II 
we establish a field-theoretic formulation for the nonlin- 



ear Langevin equation (1.1) and analyze the form of the 
propagators and vertices. Particular attention is paid to 
the symmetries of the Burgers-KPZ equation. We dis- 
cuss the time-reversion symmetry and the existence of 
fluctuation-dissipation theorems due to detailed-balance 
properties of the underlying stochastic equation of mo- 
tion. Furthermore, the Ward-Takahashi identities re- 
sulting from the Galilean invariance of the Burgers-KPZ 
equation are discussed. In Section III we proceed with 
a detailed description of our renormalization procedure 
and regularization scheme, elaborating especially on how 
we handle the d dependences of various origins appearing 
in the singular terms. We then discuss the results of our 
two-loop perturbation theory, i.e., the behavior of Wil- 
son's flow functions, the fixed points of the RG, and the 
resulting critical exponents as functions of the dimension 
d. We shall also investigate the scaling behavior at the 
dynamic roughening transition, as well as some features 
of the crossover from Gaussian to critical behavior in the 
flat phase. Finally, in Section IV we discuss our results 
and provide an outlook on future developments. In the 
Appendices A and B we present some details of the two- 
loop calculation, listing the Feynman diagrams and ana- 
lytical results for the two-point vertex functions, and the 
formulae for evaluating the integrals in the dimensional 
regularization scheme. 



II. MODEL EQUATIONS, PERTURBATION 
THEORY, AND WARD IDENTITIES 

In this section we establish a path integral formu- 
lation for the nonlinear Langevin equation Eq. (1.1), 
where particular attention is paid to the symmetries of 



the Burgers-KPZ equation. Readers who are not fa- 
miliar with these techniques, may proceed to Section 
III. A, where a formulation of the dynamic renormaliza- 
tion group is given which does not utilize the more for- 
mal methods described in this section. We note, how- 
ever, that the formulation of dynamics in terms of path 
integrals provides a quite powerful tool which allows a 
systematic analysis of higher order terms as well as gen- 
eral theorems, such as fluctuation-disspation theorems 
and Ward-Takahashi identities. 



A. The Model and Dynamic Functional 

We begin by reformulating the stochastic dynamics in 
terms of path integrals. Let us consider a coarse-grained 
interface of a d-dimensional substrate, described by a 
height function h(x, t) with x 6 lZ d , whose time evolution 
is governed by the Burgers-KPZ equation (the subscripts 
"0" will henceforth denote unrenormalized quantities) 



dh(x,t) 



An 



. ., v Q V^(x, t) + -f (V/i(x, t)Y + r?(x, t) . 

(2.1) 

The random forces f](x, t) can be taken to be Gaussian 
distributed, 



W[rf\ oc exp 



1 



4D 



d d x / dt r/(x,i) 2 \ , (2.2) 



with mean zero and short-ranged spatial and temporal 
correlations 

(r?(x,t))=0, (2.3) 
fa(x, t) V (x', 0) = 2A) 5& (x - x') S(t - t') . (2.4) 

;From the standard dynamic field theory formulation 
pl| , p2|p8|j29| ] of Langevin dynamics, upon averaging over 
the noise distribution and introducing auxiliary response 
fields h(x,t), one can turn the stochastic differential 
equation Eq.(|]l]) into a "free energy" functional (gen- 
erating functional), 



Z\j,j] 



V[h]V[ih]eaq?i J[h,h] 







j d d x j dt 


jh + jh 



(2.5) 



with the Janssen-De Dominies functional given by 
J[h, h) = J d d x J dt\ Do hh 



(2.6) 



Correlation and response functions can now be expressed 
as functional averages with weight exp |c7[/i, h]\. 
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B. Perturbation Theory 

The propagators and correlators are determined by the 
harmonic ("Gaussian") part 



Jo[h,h] = 



1 



h(k)i.Doh(-k) - [iuj + is k 2 ] h{-k) 
(*®»*®) A ®(£(j) 



(2.7) 



of the dynamic functional, with 
A(k) = 



—2D iuj + v k 2 
iuj + v^k 2 



(2.8) 



Here we have introduced k = (k, oj) and defined the 
Fourier-transformed quantities by 



h(x,t) = / / h(q,u)e 



^(qx— u)t) 



(2.9) 



using the abbreviations J ... = (2tt) d J d d q... and 
... = (27r) _1 J duj.... The interaction part J7i n t is pro- 
portional to the coupling constant Aq and given by 



Jint[h,h] 



, A 



{ki} 



{k 2 --k 3 )h(k 1 )h{k 2 )h(k 3 )6{i:k i 



(2.10) 



The corresponding diagrammatic representation of the 
response and correlation propagators and the vertex is 
shown in Fig. 1. With the above perturbation theory 
at hand we can now calculate the cumulants (...) c of 
the correlation and response functions defined by func- 
tional derivatives of F\j,j] = hi Z[j,j] with respect to 
the sources j and j, respectively: 

Gjv jW (fci; kfj] fcjy+i' •••> ^n+n) = 

= (h(ki) ... h(k^)h(kff +1 ) ... h(k$ +N )) c = 
6*+ N lnZ\j,j] 



6j(-ki)...6j(-kfi +N ) 



(2.11) 



j,j=o 



It is convenient to consider the vertex functions T 
which can be obtained from the cumulants by a Legendre 
transformation, 



T[h,h] 
where 



F\j,j}+ / d d x / dt(hj + hj) 



SF 



5F 
^3 



(2.12) 



(2.13) 



In equilibrium dynamics the linear response is defined by 
adding to the Hamiltonian an external field which couples 



linearly to the field h. For nonequilibrium random pro- 
cesses, however, the stationary probability distribution 
function (which is the analogon of the Gibbsian measure) 
is not known a priori, except for certain cases where the 
potential conditions (see below) are fulfilled and the ran- 
dom process obeys a detailed-balance condition |^o| , ^T| . 
Therefore, we define the linear response in terms of (de- 
terministic) "external" fields b added linearly to the drift 
term in the equation of motion 



V(h) -> V(h) + b. 



(2.14) 



This "external" field b corresponds to the source j in the 
generating functional. This remark should also clarify 
why h is called a "response field" (see Ref. |22]|). 

We proceed with a simple dimensional analysis to de- 
termine the engineering dimensions of the parameters 
and fields. Since the dynamic functional has to be di- 
mensionless [J] = 1, upon defining inverse length and 
time scales according to [q] — A and [oj] — vq A 2 , respec- 
tively, we arrive at the following "naive" dimensions 



[h]=K- d l^y- 3 ' 2 Dl'\ 
[h] = A-*/ 2 - 1 v- 1/2 D~ 1/2 
[A ]=A- rf / 2 +V 



-3/2 -1/2 



(2.15) 
(2.16) 
(2.17) 

Hence the engineering dimension of the effective coupling 
constant is [Ajpo/^o] = The 

space dimension at 
which the effective coupling constant becomes dimension- 
less is d c — 2. In simple cases such as a </> 4 theory this 
is the dimension which separates Gaussian from critical 
behavior. However, this simple scenario is not at all what 
one finds for the Burgers-KPZ equation. As we will see 
in the following, d c = 2 is the dimension below which the 
behavior of the interface fluctuations are described by a 
finite "strong coupling" fixed point, which is not of order 
e = d—2, but rather of order 0(1). Above d c — 2, on the 
other hand, the two-loop calculation yields an infrared 
unstable fixed point which is of order 0(e). For effective 
coupling constants less than this unstable fixed point, the 
RG flow tends to small coupling, i.e., the interface fluc- 
tuations are described by a capillary wave Hamiltonian 
Ti ex J d d x(Wh) 2 (corresponding to "Gaussian" behav- 
ior). When the coupling constant becomes larger than 
the value of the unstable fixed point, the RG flow tends to 
infinity. The corresponding strong-coupling behavior is 
hence described by a "strong-coupling" fixed point which 
is obviously not accessible by perturbative methods (of 
any finite order in perturbation theory). Resummation 
techniques such as the quite successful mode-coupling 
approach [[32| are needed. Those methods have been 
applied to the Burgers-KPZ problem in Rcfs. ||[33] p8[ . 
They allow a determination of the critical exponents and 
the scaling functions. However, little is known to date, 
whether or in what sense those self-consistency meth- 
ods constitute a controlled expansion with some possibly 
small parameter. The two-loop calculation in this paper 
provides some hints on the validity of the mode coupling 
theory in the (1 + l)-dimensional case (see Section III). 
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C. Ward— Takahashi Identities 

As usual, internal symmetries of the effective La- 
grangian, i.e. here: the Langevin equation, are important 
in determining the number of independent renormaliza- 
tion factors of the theory. The most important symmetry 
of the Burgers-KPZ equation is its behavior with respect 
to tilting of interface by an (infinitesimal) angle v (this 
corresponds to the Galilean invariance of Burgers' equa- 
tion in hydrodynamics) 



h'(x, t)= h(x + Xovt, t) + v ■ x , 
h'(x,t)= h(x + \ vt,t) . 

In Fourier space the latter equations read 
h'(q,t)=e a ° v *h(q,t)-iv-?-5(q), 



dq 



h'(q,t) 



(2.18) 
(2.19) 



(2.20) 
(2.21) 



The invariance of the generating functional of the vertex 
functions under the above transformations implies the 
following Ward-Takahashi identity 



ST 



J Jdt i\ qt ^ 
+iS(q) 



ST , . . ST ~ . . 
h(q, t) + — — -h(q,t) 



Sh(q,t) 
d ST 



dq Sh(q, t) 



Sh(q,t) 
0. (2.22) 



Functional derivatives with respect to h(q',t') and 
h(q",t"), h(q",t") yield relations between three-point 
and two-point vertex functions 



A (q'«'+qY)r k (q',t';q",t") = 
A W+q"t")T- h r{q',t';q",t") = 



d_ 
dq 



J dtT~ hJih {c{,t'-q",t"-q,t) 



q=0 



q=0 



(2.23) 



(2.24) 



or, in frequency representation 



d_ 

dui 

= r Mh(q - p> w ; q> p, o) , 

d_ 

dui 



*A (q-p)— r^(q,u/;-q,-w) = 



*A (q-p)— r^(q, w; -q, -w) = 



rEZ h (q-p,a;;-q,-«;p,0). 



(2.25) 



(2.26) 



^From the diffusive dynamics and the corresponding q 
dependence of the vertices follows the exact result (valid 
to any order in perturbation theory) 



(2.27) 



Hence the fields h a nd h do no t ren ormalize |39|. The 
Ward identies, Eqs. ( |2.25| ) and ( |2.26| ), then yield imme- 
diately that there is also no renormalization of the non- 
linearity A. Therefore we shall henceforth drop the index 
"0" for A, previously denoting the unrenormalized quan- 
tity, as we have already anticipatingly done for the fields. 



III. RENORMALIZATION GROUP THEORY 

In this section we present our results from the renor- 
malization group analysis of the Burgers-KPZ equation 
to two-loop order. We focus on the main concepts and 
results. For more details we refer the reader to the ap- 
pendices, where some intermediate results for the calcu- 
lations of the vertex functions and the renormalization 
constants are presented. 



A. General Considerations 

By now several powerful methods have evolved, and 
are frequently used, by which dynamical critical phenom- 
ena are analyzed. The method by which the Burgers- 
KPZ equation has been investigated first Q, is based 
on the dynamic renormalization group as formulated by 
Ma and Mazenko p0| . This scheme is in close analogy 
to Wilson's momentum-shell procedure 41 1, and has the 
appealing feature of being conceptually transparent and 
more "physical" than field-theoretic methods, which uti- 
lize a mapping of the sto chastic equa tions of motion onto 
a dynamical functional 



29 



nalb 
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1 1 .2 J . and concepts origi- 



■ developed for understanding quantum field theories 
]24| , [25"| . The latter schemes have the advantage of provid- 
ing a powerful tool for calculating correlation functions 
and interconnections between them in a systematic way, 
which can be crucial if one intends to go beyond one-loop 
order. 

Since the two-loop calculations will become quite com- 
plicated, and there is certainly the danger that the im- 
portant physical and conceptual points may be obscured 
by the tedious calculations, we regard it as useful at this 
point to review the dynamic renormalization group pro- 
cedure applying both Wilson's scheme and field-theoretic 
methods, to one-loop approximation. In order to keep 
our arguments as simple as possible, we also refrain from 
using the mapping onto a dynamic functional in this 
chapter, and formulate the theoretical concepts by us- 
ing jus t th e equation of motion. 

Eq. (1.1) reads in Fourier space 



h(q) = G (q)v(q) 

+ \GM J_. Vo(q+; q-)%+)%-_) . 



(3.1) 



In Eq. (3^) G (q, cj) = l/(^ q 2 — i&) denotes the "bare 
propagator", Vb(qi;q2) = — Ao qi ■ q2 is the "bare ver- 
tex", and q± = q/2 ± k, w± = uj/2 ± v. For V = 0, 
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Eq. (3.1) is just the linear diffusion equation. For Vq 7^ 0, 
the solution of Eq. ( ) may be obtained iteratively by 
a perturbation expansion in powers of Vq. For exam- 
ple, the lowest-order correction to the response function 
G = 5{h)/Sr] is 

Gi(q) = G (q) + G (q)M<l)Go(q) , (3.2) 

where Co(q) — 2Do\Go(q)\ 2 is the "bare correlator", and 

£1(9) = -^Vo(q+;q-)Go(?-)C (r+)Vo(q + ;q) (3.3) 

represents the one-loop renormalization of the "self- 
energy". Similarly, the lowest-order correction to the 
correlation function is 



C 1 (q)=2D 1 (q)\G (q) 



(3.4) 



where 

D^q) = D Q + i^ (q + ;q_)Co(g_)C (g + )K,(q + ;q_) 

(3.5) 

is the one-loop renormalization of the "noise spectrum" . 
Higher loop orders could in principle be obtained from 
iteration of the equation of motion, Eq. (3.1). But, 
since we are mainly interested in the conceptual prob- 
lems now, we restrict ourselves for the sake of simplic- 
ity to the one-loop approximation. (The two-loop con- 
tributions and their corresponding diagrammatic repre- 
sentations are presented in Appendix A and will be dis- 
cussed in Subsection B.) After performing the internal 
frequency integrals one obtains for the response and cor- 
relation functions at zero frequency w = and in the 
limit q — > 0, respectively, the following results 



d(0,0) = C (0,0) 
'" A2A, -Go(O,0) 2 



4,3 



d d k 1 
ttirYk 2 



(3.6) 



Gi(q,0)=G o (q,0) 
-G (q,0) 2 



A 2 A „ , „,,2- d 



-'I 



where only the terms of lowest order in q 2 have been re- 
tained. It is essential to note that the prefactor 2 — d in 
Eq. ( |3.7j ) is a geometry factor, stemming from the scalar 
products in the integrand, i.e., we have made use of the 
identity / k (q • k)/(fc) = (q 2 / 'd) / k k 2 f (fc) with fc =| k |. 
This d— dependent prefactor has nothing at all to do with 
the problems in perturbation theory, resulting from in- 
frared divergences, to be discussed next. Namely, the 
"naive" perturbation theory brea ks d own d ue t o the fact 
that the integral / k &~ 2 in Eqs. (3.6) and (3/7) is mani- 
festly infrared (IR) -divergent at the lower cutoff. These 
divergences are of real physical origin, and it is exactly 
those, by which one is lead to the introduction of the 
renormalization-group concept . 



In Wilson's momentum space procedure |4l| the RG 
transformation is defined as follows: 
(1) Elimination: First, one eliminates those modes within 
a small momentum shell k <E [A/6, A], where b = & > 1. 
This results in an effective (renormalized) propagator for 
the remaining modes and allows the identification of an 
effective surface tension (diffusion constant) to one-loop 
order §§ 



vq 



A 2 A (2 



A4 



KjA d - 2 l 



(3.8) 



where K4 — Sd/{2-K) d and Sd is the surface area of the 
unit sphere in d dimensions. In an analogous way one 
finds \ Wi 



A. 



Dn 



X 2 Do 

H 3 



KdA d - 2 l 



(3.9) 



Note that the expressions (3.8) and ( |3.9| ) coincide in 
d = 1, as is required by the additional fluctuation- 
dissipation theorem (FDT) valid only at this specific di- 
mension. Without loss of generality we may set the cutoff 
A to unity. 

(2) Rescaling: After this elimination-step has been per- 
formed, the resulting model has a cutoff which is reduced 
by a factor b = e > 1. In order to remove this difference 
and arrive at a form more closely resembling the origi- 
nal system, the momenta, times, and height variable are 
rescaled according to k — * ke~ l , t — > e zl t, and h — > e xl h, 
respectively. Upon requiring that the equation of mo- 
tion stays invariant under this scale transformation, the 
parameters and coupling constants transform as follows 



Do 
A 



b z - 2 v , 
b -d-2 X +z Do _ 

b* +z - 2 \ . 



(3.10) 
(3.11) 
(3.12) 



Combining these contributions results in differential re- 
cursion relations (for infinitesimal RG transformations 
where I 1) ||| 



dv w 

~dT 

dD w 
dl 

d\ w 

IT 



, o _1_ Tf ^wAv (2 - d 
z - 2 + K d 



= A, 



z - 2 X - d + K d 



4i£ V d 



4i£ 



A w [x + z - 2] 



(3.13) 
(3.14) 
(3.15) 



where the suffix "w" indicates that the running cou- 
pling constants and parameters are obtained within the 
framework of Wilson's momentum shell method. (We 
shall later introduce another set of flow functions in the 
context of renormalized field theory.) Note that there 
is no perturbative correction to the flow of A w , which 
can be understood as a dir e ct con sequence of the Ward- 
Takahashi identities Q2 23| )-( |2^ ). 



G 



(3) Fixed point and critical exponents: The last step is 
to determine the exponents z and x such that the pa- 
rameters v w and D w are not changed, for we have then 
obviously arrived at a scale-invariant situation. We de- 
fine an effective coupling constant by 



g w = Kd 



\ 2 n 

W w 

4l/ 3 



(3.16) 



and obtain for the corresponding flow: 
dg w (l) 



dl 



(2 - d)g w (l) 



^4(0 



(3.17) 



^From these recursion relations one can solve for the 
RG fixed points, i.e., points in parameter space which 
are invariant under sc ale tr ansformations. Thus setting 
dg w (l)/dl = in Eq. (3.17) yields the "Gaussian" fixed 
point (?q = 0, and also one non-trivial fixed point 



fJi 



d(d-2) 
2(2d— 3) ' 



(3.18) 



which is positive for < d < 3/2, diverges at d = 3/2, 
then rises again from minus infinity, crossing the zero 
axis at d = 2, and becoming positive for d > 2 (Fig. 2). 
In the physical region, g > 0, g\ turns out to be sta- 
ble ("attractive") for d < 3/2, thus providing non-trivial 
(strong-coupling) scaling behavior, while for d > 2 the 
Gaussian fixed point is stable. Hence, depending on 
its initial value, g w {l) will either flow to <?q, or to the 
strong-coupling "fixed point" g w — > oo. Thus for d > 2 
a dynamic phase transition ("roughening transition") is 
found, governed by the unstable fixed point g\. The fact 
that there is no finite positive fixed point in the range 
3/2 < d < 2 will turn out to be an artefact of the one- 
loop approximation (see Subsection B). 

Upon setting dv^jdl = and dD w /dl = 0, and g w — 
g*, one finds zq = 2 and xo = (2 — c?)/2 at the Gaussian 
fixed point, while at any non-trivial fixed point g* ^ 
the critical exponents become 



z = 2 



X 



d-2 



d-2 



(3.19) 
(3.20) 



Note that Eq. (|3T5]forces z + x = 2, if A ^ 0. Thus the 
exponent identity (1.6) holds for any finite fixed point, 
which implies that knowledge of either the flow of v w or 
D w suffices in order to determine the critical exponents 
in its vicinity. 

(As already mentioned, the flow obtained from these 
one-loop equations still contains some artefacts, which 
are due to the low order of perturbation theory consid- 
ered. Yet, our concern in this chapter is rather a com- 
parison of the various methods available to determine the 



RG flow of the parameters and couplings than a detailed 
description of the physical behavior.) 

Let us now compare these results with the various field 
theoretic techniques using different regularization proce- 
dures. In field theory, the cut-off A is set to infinity, thus 
leading to ultraviolet (UV) divergences above a certain 
critical dimension d c (in our case, d c = 2). There are ba- 
sically two regularization methods, by which these inte- 
grals are assigned meaningful values, that are frequently 
used in field theories: (i) reintroducing the cut-off, and 
(ii) dimensional regularization. In the former case the 
original UV singularities appear as logarithms or pow- 
ers of A, in the latter as l/(d c — d) poles. The criteria 
of choosing either of them are mostly guided by conve- 
nience, i.e., the amount of effort by which the correspond- 
ing integrals may be evaluated. Using those more formal 
field-theoretic methods one has to be careful, however, 
especially if one deals with massless theories, as is the 
case here! Let us explain why some specific attention is 
required (see also Ref. [ p5[ , p. 161). In dimensional regu- 
larization |E3,E4| one has 



d d p 



0. 



(3.21) 



which is a (non-trivial for even m) consequence of the 
dilatation property of the integrals defined here by di- 
mensional continuation. This result can be regarded 
as a cancellation between infrared (IR) and ultraviolet 
(UV) divergences. This often convenient property of inte- 
grals within the dimensional regularization scheme may, 
however, have quite dangerous consequences: In a field 
theory involving massless fields (for which the Burgers- 
KPZ equation is an important example), the latter gen- 
erate IR singularities, which, again, have the signature 
of l/(d c — d) poles! Thus, contrary to calculations using 
a finite cut-off, the dimensional regularization scheme is 
in d anger of mixing UV and IR poles, as is obvious from 
Eq. (3.21). It is, however, crucial to clearly separate the 
UV divergences, and avoid any mixing with the IR sin- 
gularities, as will be explained below. This really may 
be the source of many mistakes, if this method is ap- 
plied without sufficient precaution. Thus it is essential 
to employ an IR cut-off in the integrals. In the case of 
the Burgers-KPZ equation, we cannot simply introduce 
a mass term, however, because this would immediately 
violate the Galilean invariance discussed in Section II. C. 
Therefore we have to evaluate all quantities explicitly at 
finite external wave vectors or frequencies. For the cal- 
culations in Appendix B, we have chosen the convenient 
normalization point NP: q = 0, ioj/2v = ff 2 . 

In order to make our arguments clear we want to 
very briefly review some of the fundamental ideas behind 
renormalized field theory and the connection to critical 
phenomena |2jj2^,^]. From the dimensional engineer- 
ing in Section II and simple power counting arguments, 
it follows that perturbation theory fails to describe the 
critical theory below the critical dimension d c because of 
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IR divergences (see above). Above d c , on the other hand, 
the integrals become UV divergent. These UV singular- 
ities may be handled very effectively by using the well- 
established technology of quantum field theory, namely 
by introducing renormalization constants into which the 
UV poles are absorbed. From these the Callan-Symanzik 
equations can be inferred, which describe the dependence 
of the renormalized quantities on the renormalization 
scale. Already by a heuristic argument (see Ref. p2[ , 
p. 271), the intimate relation of the IR and UV singular- 
ities may be explained as follows. As mentioned above, 
we shall have to carefully employ a finite IR cut-off, say 
a mass parameter fj, (which in our case is rather the ex- 
ternal frequency). The physical IR divergences in which 
we are primarily interested, manifest themselves in the 
domain q/ fj, 3> 1, and q/A <C 1, where q are the ex- 
ternal momenta. By construction the latter condition is 
automatically satisfied in renormalized field theory. The 
former can also be read as q — > oo for fixed [i. Note 
that such a scaling argument is possible, only when the 
loop-integrals are arranged in such a way that they are 
insensitive to their upper limits (this is what one does in 
renormalizing the theory). More specifically (and tech- 
nically), a change in the normalization point fi — > bfi 
of the theory may be interpreted as a scale transforma- 
tion q — > q/b, and solving the Callan-Symanzik equation 
provides the corresponding RG flow equations. Again 
the fixed points with respect to these (infinitesimal) scale 
transformations are investigated, and the fixed point val- 
ues of Wilson's functions, the so-called "anomalous di- 
mensions", yield the critical exponents. 

Let us explain this for the case of the Burgers-KPZ 
equation. In order to renormalize the theory (i.e., re- 
move the UV divergences) one introduces renormalized 
parameters, which are related to the "bare" parameters 
through renormalization factors containing all the UV 
poles. Before turning to the one-loop results, some gen- 
eral features of the renormalization, which are valid to 
every order in perturbation theory, should be noted. As 
discussed in th e pre vi ous s ection, the Ward- Takah ashi 
identities Eqs. fl2.23| )- (|2.26|) together with Eq. (|2.27|) en- 
sure that the renormalization involves only two indepen- 
dent renormalization factors, namely, the renormaliza- 
tion of the noise amplitude Dq and the surface tension 
Uq. Hence we define renormalized parameters by 



D = Z D D 



V = Z v Vq, 



(3.22) 



and determine these by the following normalization con- 
ditions for the singular parts of the two-point vertex 
functions 



^ t-> / \ i sine 

^2 r ^(q^)l NP - 



-2D. 



(3.23) 
(3.24) 



The normalization point (NP) is conveniently chosen at 
q = and iiojlv = [i 2 , and is an arbitrary length 



scale. We remark that this is not a minimal-subtraction 
prescription, as was applied by Sun and Plischke [p7| , 
where just the residues of the l/(d c — d) poles (in the 
dimensional regularization scheme) would be included in 
the Z factors. Rather we retain the complete dependence 
on the dimension d in the geometric prefactors originat- 
ing in the an gular i nteg rations [see the discussion follow- 
ing Eqs. (0) and @]. 

The corresponding Wilson functions , (d and the /3 
function (3 g (see below) for the effective coupling constant 



,9() 



AoA) 



(3.25) 



permit us to study the renormalization-group flows of the 
renormalized parameters and coupling constant D, v and 
g. Note again that due to the Ward-Takahashi identities 
there is no renormalization for the amplitude A of the 
three-point vertex to all orders of go, i.e., Z\ = 1. Hence 
one gets for the renormalized effective coupling constant 



g = Z g9o C d fi d - 2 with Z g = Z D Z- 3 , (3.26) 

where we have absorbed the geometry factor C d = 
K d T(d/2)T(2 - d/2) = T(2 - d/2)/2 d - 1 n d / 2 in the defi- 
nition of the renormalized coupling. 

To one-loop order the bare vertex functions read 



Re 



hh 

= -2D 
^%(q, W ) sl "» | q =o= 



iu /2vq 



X^Dq d-2 



Re 



iuj/2i>Q + k 2 



(3.27) 



(3.28) 



where the <i-dependent prefactor in the latter equ ation 
is of exactly the same geometric origin as in Eq. (3.7) 
above. In d = 1 the results for the amplitude of the noise 
spectrum and the surface tension again become identi- 
cal as is required by the fluctuation-dissipation theorem 
valid in d = 1. This important feature would have been 
lost in a minimal subtraction prescription! j4^| . Apply- 
ing dimensional regularization, the above integrals yield 
at the normalization point iuj/2v = [i 



2. 



d-2 



k M 



k 2 



d-2 



(3.29) 



The l/(d — 2) = 1/e pole corresponds to a In A in the 
cutoff-regularization scheme, and to the momentum shell 
integral oc / in Wilson's scheme. At this point we empha- 
size that this 1/e pole has to be subtracted in order to 
renormalize the theory (i.e., remove the UV divergences 
corresponding to In A at the critical dimension d c = 2). 
Let us now demonstrate once more why the dimensional 
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regularization method (with minimal subtraction) is a 
dangerous procedure unless carried out with considerable 
precaution. Namely, one might argue that the 1/e pole 
for the noise amplitude is cancelled by the prefactor d— 2 
and hence does not have to be incorporated in the renor- 
malization factors |27|| . In our opinion, that would be 
incorrect! This procedure would definitely leave UV di- 
vergences in the theory, as becomes obvious in the cut-off 
regularization scheme, where In A terms would survive. 
In summary, one should keep in mind that the 1/e poles 
in the dimensional regularization scheme are essentially 
nothing else but a quite convenient way of keeping track 
of the UV (In A) divergences at the critical dimension 
d c = 2 (the simplification is that the resulting integrals 
are much easier to carry out in the dimensional regular- 
ization method); one should thus strictly avoid mixing 
these "artificial" dimensional factors with others origi- 
nating from purely geometrical properties of the integrals. 
Again, this might constitute a troublesome trap unless 
sufficient attention is paid to this issue. Henceforth, and 
in the Appendix B, we try to emphasize this by clearly 
distinguishing between "d" and "e" . 



Using the definitions (3.22), the normalization condi- 
tions (3.23), (3.24), and the one-loop results (3.27) and 



(3.28) for the singular parts of the two-point vertex func- 
tions, one arrives at 



Zn = l- 



Z v = 1 



goCd^ 



d-2\ g C d pL e 



Upon defining Wilson's flow functions by 



t 9 

6[l 



D : 



r d 



In Z u 



and the [3 function 
d 



Pg = A* 



On 



g(d-2 + C D - 3( v ) 



we find the following the one-loop results 
( D = -g + 0(g 2 ), 



d-2 



9 + 0(g 2 ), 



, , 2(2d-3) 



(3.30) 
(3.31) 

(3.32) 

(3.33) 

(3.34) 
(3.35) 

(3.36) 



Note that these "physical" quantities do not contain the 
artificial e factors any more! 

Running parameters and coupling "constants" are now 
defined by (see Subsection B) 



l^ = ( D (l)D(l) 



(3.37) 
(3.38) 
(3.39) 



We remark that these flow equations are related to the 
flow equations in Wilson's scheme by the following re- 
placements [see Eqs. ( 3.13 ), ( 3.14 ), and ( 3.19| )]: 



ln(l//H« : 



!/-> *Av = ve^ 1 . 



(3.40) 
(3.41) 
(3.42) 



By determining the zeros of the /3 function ( |3.33|) we find 
the same fixed points g$ and g\ as in Wilson's scheme 
[Eq. (3.18), Fig. 2]; note, however, that we have incor- 
porated somewhat different geometry fact ors in the def- 
initions of the renormalized coupling g (3.26) and g w 
fl3.17| ). As will be demonstrated in the next subsec- 
tion, the critical exponents are given by z — 2 + £* and 
X = 1— (rf+C|) — C)/2j respectively. This yields the iden- 
tical results as the momentum shell procedure, demon- 
strating (to this order) that both schemes are equivalent 
methods in order to find the universal critical behavior. 

A very important final remark is in place here. In ei- 
ther scheme, the above one-loop results were obtained by 
evaluating the integrals at fixed dimension. If an expan- 
sion near d = 2 were applied, hardly any answer would 
have been found for d < 2, due to the divergence of the 
fixed point gl at d = 3/2. This is precisely the regime 
which is addressed by Sun and Plischke p7fl in their re- 
cent two-loop calculation. They do indeed find another 
finite fixed point at d = 2 with their approach, but it 
is difficult to see how this result may be consistent with 
the implicit assumption of the e expansion that any non- 
trivial fixed point be of order e. Our own findings within 
the two-loop approximation, as explained below, rather 
seem to indicate that a "full" 2 — e expansion cannot 
be performed consistently, in accord with the trend al- 
ready seen on the one-loop level. The situation is en- 
tirely different for d > 2. Here one may follow the ideas 
exploited in the study of the 0(n) nonlinear a model 
p6[ , and perform an e expansion above the critical di- 
mension d c = 2. "Naive" power counting suggests that 
the theory should not be renormalizable for d > d c , which 
would mean that the critical behavior could no more be 
inferred from studying the UV limit of the theory, which 
would be ill-defined. However, there appears a new fixed 
point of order e, which is IR-unstable, i.e., W-stable, 
meaning that it governs the large-momentum behavior. 
Hence the theory is renormalizable despite the "naive" 
power-counting arguments. This is precisely what allows 
for a description of the critical properties of the model, 
for the IR-unstable fixed point physically corresponds to 
a second-order phase transition in the system. Indeed, 
if one inse rts the one- loop fixed point value ( |3.2C| ) into 
Eqs. ( 3.21 ) and ( [3.22] ) for the critical exponents at the 
roughening transition, one finds that z c = 2 + 0(e 3 ) and 
Xc = + C(e 3 ), in accord with the scaling argument by 
Doty and Kosterlitz @. 
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B. Two—loop Results 

In this section we return to the fi eld- theoretic method, 
based on the dynamic functional ( |2.6|) , and present our 
results from the two-loop calculation. We shall primar- 
ily focus on the analysis and discussion of the results. 
For some details on the calculations and technicalities 
we refer the reader to the Appendices. There is, how- 
ever, one important point that we have to mention here. 
Namely, to two-loop order the situation becomes addi- 
tionally complicated with respect to the previous discus- 
sion, due to the fact that our massless field theory pro- 
duces new IR singularities at d — 1 and d — 3 (compare 
Appendix B). These can only be handled by a "partial" 
e expansion; i.e., we keep all the geometry factors, as 
explained at length above, and expand only in those d— 
dependent factors stemming from the integrals as calcu- 
lated in the dimensional regularization scheme, indicated 
with "e" instead of "d" in the Appendix B. This is, we 
admit, a most subtle procedure, and together with the 
already mentioned features, it marks the essential differ- 
ence to the approach by Sun and Plischke ^7|. To our 
opinion, however, this is the only possible method in or- 
der to approach the regime d < 2 consistently, keeping 
the exactly known properties at d — 1. For d > 2, we 
shall eventually perform a "full" e expansion in the fi- 
nal results, where all the problems cited in the previous 
subsection have already been accounted for. 

^From the two-loop expressions f or th e s ingul ar parts 
of the two-point vertex functions (5.23), ( 5.24 ) i n Ap - 
pendix A, the normalization conditions ( 3.22 )— ( 3.24 ), 
and the integrals in Appendix B, we find the following 
results for the renormalization constants 



d-2, 



d 2e 



+(d-l)2£+0<$), 



(3.43) 



Z v = 1- 



d-2 



9 ^ + (d-l) 9 f 
e 2e 



d - 2 g 2 a 
d 2e 

d - 1 9q 



F„(d)+O(g 3 ), 



(3.44) 



d 16e 

where go = ggCd^, and we have defined 

F v {d)= 4 - 2(6 - d)/ 00 (2) + 2rfJ 10 (2) 

+21/ n (2) - J ai (2) - 7J 22 (2) + 2/ 32 (2) 

+4M2)-4A 2 (2)-^ln^4 (3.45) 

The parameter integrals I rs (d) and I rs {d) are defined in 
Appe ndix B. The corresponding Wilson flow functions, 
Eqs. ( 3.32 ), are given by 



C D = - g - (d - 1) !. 2d - 1 V + Oto 3 ), (3-46) 



d-2 (d-l)(d-2) 2 
Cu = —1—9 H 1 9 



d- 1 

8d 



g 2 F„(d) + 0(g 3 ) , 



(3.47) 



/ 2(2d-3) (d-l)(5d-7) _ 2 
(3 g =g\d-2 g g 



d-l 
8d ' 



2 F u (d)+0(g 3 ) 



(3.48) 



(Note that there are no e-dependent terms left in these 
expressions, specifically, there are no 1/e poles, which 
constitutes a very non-trivial check to the calculations, 
along with the fact that at d = 1 the FDT is fulfilled, as 
required, see Appendix A.) 

Most interestingly, the two-loop contributions to these 
£ functions vanish at d = 1! That is, there are no sin- 
gular contributions to the two-point vertex functions in 
one dimension] This is clearly a very valuable fact for 
the purpose of justifying a self-consistent approximation 
as the mode-coupling approach, where vertex corrections 
are neglected |54j . In terms of the fixed points, this means 
that at d = 1 the strong-coupling fixed point in the two- 
loop approximation is unaltered with respect to the one- 
loop result: g\ — 1/2, and, of course, the critical expo- 
nents are not modified either, z = 3/2 and \ = 1/2- This 
is reassuring, for these values already follow from a com- 
bination of t he FDT (t hus C * = ( p) a nd the exponent 
identity (Q, see Eqs §M ) an d j|j5f ) below. 

The flow functions ( 3.46 ) and ( |3.47|) permit us to study 
the renormalization-group flow f or th e renormalized ef- 
fective coupling constant g, Eq. ( 3.33| ), as a function of 
the dimensionality of the growth problem. The general 
features of this flow, and the behavior of the ensuing three 
fixed points <?q , g\ and g\ as functions of the dimension 
d are (Fig. 3): 

(i) Just below the borderline dimension d c = 2 there are 
two (non-negative) fixed points, the Gaussian fixed point 
<7o = 0, which is unstable, and one strong coupling fixed 
point g*{. The "weak-coupling" fixed point gg = 0, de- 
scribing a smooth interface, is IR-unstable. Hence for 
d < 2 the RG flow always tends to the strong-coupling 
fixed point g\, describing a rough surface. (Below d = 1 
there is an additional fixed point in the physical region, 
whose value diverges for d — > 1. This would constitute 
an unstable fixed point in the flow of the coupling; it ap- 
pears, however, rather doubtful to extrapolate the results 
of our two-loop calculations beyond d = 1, as we had to 
apply the above-mentioned "partial" e expansion.) 

(ii) At the critical dimension d c = 2 the strong-coupling 
fixed point g\, as obtained from the two-loop calculation, 
tends to infinity. It is not clear, whether this divergence 
of the strong-coupling fixed point at the critical dimen- 
sion d c = 2 is just an artefact of the two-loop calcula- 
tion or a general feature of any finite order in perturba- 
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tion theory. We suppose that it indicates that no finite 
strong-coupling fixed point will emerge to any order in 
the perturbation expansion. In this sense, there is non- 
perturbative strong-coupling behavior for alld > 2. There 
is still a considerable amount of controversy about the 
existence of an upper critical dimension (not to be con- 
fused with the critical dimension d c = 2 obtained from 
power counting arguments), at which mean field values 
are recovered. Actually, in view of the analogy to the 
nonlinear a model, the critical dimension d c — 2 can be 
regarded as the lower critical dimension, since for d c < 2 
no dynamic roughening transition takes place [see (i)]. 
The random directed path on a Cayley tree [Q repre- 
sents a possible candidate for a high-dimensional limit. 
The existence of a finite upper critical dimension is sup- 
ported by an expansion by Derrida and Cook |45|] , who 
insert finite sections of the high-dimensional lattices in 
place of the nodes of the Cayley tree. The validity of 
such an approach has been criticized by Fisher and Huse 
Ji2| |, however. Also, numerical simulations [|l4|-[l6| sug- 
gest that there is no upper critical dimension. From our 
two-loop results it is not possible to draw any decisive 
conclusion about the existence of an upper critical dimen- 
sion. Even higher loop orders are most likely not to be 
very useful for deciding upon the question of the existence 
of an upper critical dimension. What is really needed is 
a systematic expansion in 1/d, or some other controlled 
expansion capable of taking into account an infinite set 
of diagrams. A very recent mode-coupling analysis fl38[ , 
which makes no assumptions about the line shape, but 
is still based on an uncontrolled approximation, indeed 
yields z < 2 for all d. 

(iii) Above d c — 2 our two-loop calculations support the 
existence of three fixed points. There are two IR sta- 
ble fixed points, whose domain of attraction is separated 
by a critical fixed point g% — g^, which is IR-unstable 
(hence UV-stable). The critical fixed point g* describes 
a dynamic phase transitions and is accessible by pertur- 
bation theory (it is of order e = d — 2). For g < g* 
the RG flow tends towards the weak-coupling fixed point 
5wc = 9a = 0, describing a smooth interface. For effec- 
tive coupling constants larger than g* the flow leads to 
an IR-stable strong-coupling fixed point g\ = g^ c - This 
strong-coupling fixed point seems to be inaccessible by 
a perturbational approach. To two-loop order one thus 
finds g* c = oo. 

This scenario is similar to the results obtained for the 
0(n) nonlinear a model. In a 2 + e expansion pq 1 ] one 
finds a non-trivial zero of the corresponding [3 function. 
This IR-unstable fixed point defines the critical temper- 
ature, in exactly the same way as the above IR-unstable 
fixed point in the roughening problem defines the criti- 
cal coupling, at which a dynamic phase transition from 
a smooth to a rough surface takes place. This analogy 
becomes even more apparent if one considers the map- 
ping of the Burgers-KPZ equation onto the statistical 
mechanics of directed polymers in random media. For 
d > 2 Imbrie and Spencer M have shown rigorously 



that the polymer undergoes a continuous transition from 
a low temperature pinned phase to a high temperature 
phase where the disorder is irrelevant. The above found 
critical fixed point controls this transition. 

We now proceed to a discussion of the behavior near 
the different fixed points for d > 2. In order to relate 
the values of the £ functions at the fixed point to the 
dynamic exponent z, and the roughness exponent \i it is 
convenient to consider the two-point correlation function 
Ch/i(q, lo), which acquires the following scaling form 



C hh (q,u;) = q- d -^- z C(u;/q z ), 
or equivalently (see Section I) 

C hh (x,t) =x 2 *C(t/x z ). 



(3.49) 



(3.50) 



We want to analyze how this scaling form and the expo- 
nents for the correlation function are related to the re- 
sults obtained from renormalized field theory. Since the 
bare vertex functions are independent of the arbitrary 
momentum scale [L introduced in the RG procedure, one 
finds the following RG equation for the two-point vertex 
functions T(fi, v, D, g, q, uj) 



d_ 







r(.)=0, (3.51) 



where we have introduced the abbreviation (.) = 
(p., v, D, <?, q, lo). The RG equation is readily solved with 
the metho d of characteristics. The characteristics a(l) of 
Eq. (3.51) define running coupling "constants" and pa- 
rameters into which these transform when /i — > /i(Z) = fjd. 
They are given by the solutions to the flow equation of 
the coupling, Idg/dl = P g (l), and the first-order differen- 
tial equations for the parameters a — D,i> 



(3.52) 



with the initial conditions D(l = 1) = D and v(l = 1) 
v, namely 



a(l) = aexp 



dp 



Up) 



(3.53) 



Applying a dimensional analysis one finds that and 
1TV, have dimensions D and Ufi , respectively. Hence the 
solutions of the RG equations read 



9(1) 



r^(-) = £>(0r* 



hh 



fiV p?l 2 v(l) 



,9(1) 



(3.54) 
(3.55) 



Since the two-point correlation function is related to the 
two-point vertex function by C = T 



', one gets 

at a fixed point, with the matching condition q/ fd = 1, 



hh' I hh I 
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2+c; 



(3.56) 



Hence we arrive at the following, already mentioned re- 
lations 



X = l 



c - (h 



(3.57) 
(3.58) 



z = 2 + C- 

At the weak-coupling fixed point, <?q = 0, we have 
X = 1 — d/2 and z = 2. At any non-zero fixed point 
g* =^ one gets 



d-2 + &-3C = 



(3.59) 



from f3 g (g*) = 0. Note that this exponent identity results 
ultimately from the Galilean invariance of the Burgers- 
KPZ equation. Hence, we find 



X = ~W ! 

2 = 2 + C 



(3.60) 
(3.61) 



^From these relations one can easily infer the exponent 
identity x + z = 2 already mentioned in Section I. 

Let us now investigate the two-loop results at the IR- 
unstable (critical) fixed point g*. One can show that 
F v (d) = 8 + 0(e), where the 0(e) coeffficient can also 
be d eterm ined numerically F v (d) « 8 — 4.0797e [see 
Eq. ( 5.38j )l . Hence, one finds that in a consistent ex- 
pansion in both g and e 



Cd = -9~ -9 2 + 0( 5 3 ,5 2 e, ff e 2 ,e 3 ), 



e 1 
2 9 ~ T 



C, = x5-x.9 2 + 0(5 3 ,5 2 e,5e 2 ,e 3 ) 



(3.62) 
(3.63) 



where e = d - 2. With Eq. fl3.33| ) this yields for the /3 
function 



— .. I . ^ + 3e _ | /n/ 3 ^ 2 



.9 



:.g + 0(^,^e, ff e^,e J ) 



(3.64) 



which is identical to the one-loop function since all the 
0(g 2 ) corrections cancel! The resulting critical fixed 
point is g* = e — 3e 2 /2 + 0(e 3 ), and there is no finite 
strong-coupling fixed point. Inserting this result for the 
unstable fixed point into the C functions one realizes that 
all 0(e 2 ) corrections cancel, i.e., 



G = -e + 0(e 3 ). 
C = + 0(e 3 ) . 



(3.65) 
(3.66) 



Hence our two-loop calculations confirm the results by 
Doty and Kosterlitz Jig] , that z c = 2 and \c — at 
the roughening transition. Note that in our perturbation 
expansion this result is due to a most remarkable and not 
at all obvious cancellation of very different contributions. 
This is a very reassuring feature of our method again, at 
least if one follows the considerations in Ref. lfl8| . 



Now let us reverse the argument, and assume that z c = 
2 is a n exa ct res ult at the critical fixed point. Then, with 
Eqs. (3.59) and (3.61) one gets at the weak-coupling fixed 
point 



(d-2) + ( D (g* c ) =0, 



(3.67) 
(3.68) 



Therefore, it should also be possible to determine the 
value of the critical fixed point from one of the latter 
equations. In fact, one obtains from Eq. (3.67) the same 
critical fixed point as from the zero of the /3 function. 
However, in order to determin e the value of g* up to 
terms of order 0(e 3 ) from Eq. ( |3.68D one needs the co- 
efficient of the 0(g 2 ) term and the 0(g 3 ) term of £„ up 
to terms of order 0(e 2 ) and 0(e 1 ), respectively. For the 
calculation of the fixed point value from the zeros of the 
(3 function, we had to know the Q(g 2 ) term of Q v up to 
terms of order 0(e), only. From Eq. ( 3.47 ) one obtains 
upon including an unknown three-loop term 



, 4.0797 
+ — *).</■ 



(3.69) 



where A is a constant of order 0(1). Eq. (3.68) yields 
for the critical fixed point: g* = e — 3e 2 /2, only if 
A = —(3/2 + 4.0797/8), i.e., one can determine the 
three-loop correction for the renormalization of the sur- 
face tension v. 

Let us comment on the d dependence of the critical 
fixed point. Keeping the full dependence of the C fu nction 
for the noise amplitude on d, one obtains from Eqs. (3.46) 
and (3.67) for the critical fixed point 



-1 + y/l + 4(rf - l)(2d -l){d- 2)1 d 
2(d- \){2d- l)/d 



(3.70) 



which in the limit of large d reduces to a finite value 
g* — > l/v2. It would be interesting to sec how this d 
dependence of the fixed point compares with numerical 
simulations. Note, however, that the scale of g* is non- 
universal and depends on the precise definition of the 
coupling g. 

We close this section by a discussion of the RG equa- 
tions for d > 2, and its implications on the kinetic rough- 
ening transition described by the IR-unstable (repul- 
sive) fixed point g*. We remark that a detailed analysis 
of the crossover scaling behavior in the weak-coupling 
regime at the one-loop level has been given by Natter- 
mann and Tang j2(J. We have shown above that the 
exponents at the transition are given by z c = 2 + 0(e 3 ) 
and Xc = + 0(e 3 ). In the smooth phase (i.e., for cou- 
pling constants g < g*) the effective coupling approaches 
zero. Therefore, the long-distance and long-time prop- 
erties in this phase may be obtained from a perturbative 
RG study. (In the rough phase a perturbative expansion 
is not possible due to the strong-coupling behavior.) 

The proximity to the critical point can be measured in 
terms of the control parameter S = {g* c ~9)lgt^ an d, quite 
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analogously to the treatment of the nonlinear a model at 
its critical point p6fl , we define a correlation length £(g) 
via the solution of the differential equation 



(3.71) 



Since £ has the dimension of a length this can also be 
written as 



9 a is d 



^) = M-y /(d - 2) 



exp 



£(s,m) = o, 



dg' 



(3.72) 



1 



/W) {d-2)g> 

(3.73) 



Close to the critical fixed point we can write 



dg p g (g(l)) [g(l) - g* c ]P' g (g* c 

the solution of which is given by 



(3.74) 



(3.75) 



where (3'(g*) = d0 g (g)/dg\ g *. Hence the critical expo- 
nent for the correlation length is determined by the first 
derivative of the function at the critical fixed point 
v = —l//3' g (g*). Since we have shown above that all 
two-loop contributions to this function cancel, we find 
that the correlation length exponent is 



l/is = e + 0(e 3 ) 



To one-loop order, this result was already obtained by 
Nattermann and Tang (20), who also discuss how the 
above scaling is modified to a logarithmic f orm a t d c = 2. 

Let us now solve the RG equations, Eq. (|3.5l| ), in such 
a way that the scaling behavior at the dynamic phase 
transition becomes apparent (see also Ref. |2C|l ). Dimen- 
sional analysis tells us that 



^hh^,D,g,q,u) = DTtji -,-3- 



r w,0*,f>Aff,q.f) = [i 2 v^~ Ylh [-,,-57, 



q UJ 



(3.77) 
(3.78) 



Thus the solution of the RG equation can be written as 

^ 2 (g) 



r^(-) = D( 9 )r h - h <?£(<?), 



r- hh (-) = ag)- 2 Hg)r hh <?£(<?), 



via) 



(3.79) 
(3.80) 



with 



v{g)=v ^{-Slm U9 ' 



(3.81) 
(3.82) 



In summary, one obtains the the following dynamic scal- 
ing behavior of the correlation function 



C(<y,q,«)= M(£( ff )K 4 ( 5 )(7 ( q t(g), 



(3.83) 



where we have defined the quantity 



M(((g)) = exp 



dg 

W) 



{( D (g)~2Ug)}> ■ (3-84) 



We have also introduced a characteristic frequency by 



"> c (£(ff)) = A* £ ex P1 " 



dg 

W) 



Ug) 



(3.85) 



Close to the critical fixed point the latter equation re- 
duces to ui c (£(g)) oc £ _Zc , with the dynamical critical 
exponent z c = 2 + Q. One should note that this dy- 
namic scaling form is quite analogous to the correspond- 
ing scaling in the nonlinear a model [ p6| , and the same 
conclusions drawn there apply for the roughening transi- 
tion as well. Let us contrast the above scaling behavior 
with the well-known scaling behavior of a standard 4> 
model close to its critical point. The 4> A model depends 
on two coupling constants, the coefficient of the </> 2 term 
which plays the role of the control parameter (reduced 
temperature) corresponding to the reduced effective cou- 
(3.76) phng S defined above, and the coefficient of the (j) 4 term 



which has no equivalent here. The scaling form of the 
4> model is obtained at the IR-stable fixed point, e.g. 
obtained by an e expansion near the upper critical di- 
mension 4. In the present case, however, there is not 
only scaling behavior at the critical fixed point, but in 
the entire "smooth" phase. This behavior, Eq. (3.83), 



is most similar to that of 0(n) symmetric models in the 
low temperature phase, and more closely resembles the 
form known from crossover scaling behavior J20| ]. In or- 
der to study the crossover scaling behavior in the "weak 
coupling" phase, we use the following "resummed" ex- 
pressions for the C and functions 



Cd 



g 



1 - 3/2 ff ' 



> _ g_ ( g 

^ " d \ 1 - 3/2g 



9 M 



2 + 3e 



(3.86) 
(3.87) 
(3.88) 



which are identical to the expressions above, but also 
take into account that the fixed point values can be 
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equally w ell ob taine d fro m the zero of the (3 function 
and Eqs. (^67|) and flUf ). 

First, we consider the crossover of the correlation 
length as a function of the effective coupling g. As de- 
picted in Fig. 4, the correlation length crosses over from 
£ oc g 1 ^ at small couplings to £ oc \g — g*\~ u as g ap- 
proaches the critical coupling g* . The crossover function 
shown in Fig. 4 is universal, i.e. all the non-universal 
scales may be absorbed in an amplitude £o for the corre- 
lation length. Specifically, the location of the crossover 
Scross can be obtained from Fig. 4. Similar crossover be- 
havior is found for the noise amplitude D and the surface 
tension v, as shown in Fig. 5. The asymptotic behavior 
at the critical fixed point is given by 



D(g)cx\g-g* c \-' 
%)°c|<7- 9c T. 



\9-9* c \~ 



(3.89) 
(3.90) 



The noise amplitude shows a crossover from an expo- 
nential increase D{g) oc e 9 ^ at small values of g to 
a power law divergence D{g) cx \g — g*\~ eu as g ap- 
proaches the critical coupling g*. The renormalized sur- 
face tension v(g) crosses over from exponentially decreas- 
ing v(g) cx e~ 9 / 2 to a constant at criticality. Finally, the 
crossover behavior of the characteristic frequency is de- 
picted in Fig. 6. It crosses over from u) c (g) oc g~ 2 ^ € e~ 9 ^ 2 
at small couplings to oj c {g) cx \g — g*\ as g approaches 
the critical coupling g* . 

It would be interesting to test the validity of the above 
crossover scaling, including the from of the crossover 
functions, by numerical simulation. 



IV. SUMMARY AND CONCLUSIONS 

In this paper we have given a systematic analysis of the 
Burgers-KPZ equation in d + 1 dimensions by dynamic 
renormalization group theory. We have paid special at- 
tention to the interconnections of the various renormal- 
ization group techniques. 

Let us briefly summarize our main conclusions. For the 
roughening dynamics as described by the Burgers-KPZ 
equation we find the following results: (1) The roughen- 
ing transition of the Burgers-KPZ equation is understood 
in terms of an IR-unstable fixed point above the critical 
dimension d c = 2. The value of this critical fixed point 
is accessible by a 2 + e expansion, similar to an analo- 
gous expansion for the nonlinear er model. The critical 
properties of this transition are characterized by one in- 
dependent exponent. Doty and Kosterlitz have argued 
on the basis of a standard scaling argument that the dy- 
namic exponent z c at the dynamic roughening transition 
should be exactly equal to 2. From our two-loop cal- 
culation we find that z c = 2 + 0(e 3 ), which supports 
their considerations. We have further analyzed the scal- 
ing behavior near the transition and in the smooth phase. 
Introducing a control parameter S — (g* — g)/g* 1 which 



measures the distance to the critical point, one can in- 
troduce a correlation length £ and calculate the corre- 
sponding exponent v. We find that v = 1/e + 0(e 3 ). (2) 
Below the borderline dimension d c — 2 there appear two 
fixed points, an IR-unstable Gaussian fixed point and 
an IR-stable "strong-coupling" fixed point. Note that 
the strong-coupling fixed point is not of order 0(e) and 
there is apparently also no other known small parame- 
ter. Hence, the phase described by this strong-coupling 
fixed point is not accessible by a controlled perturba- 
tion theory, and one needs non-perturbative methods 
(e.g. mode-coupling theories) to access the physics of 
the rough phase. The "weak-coupling" fixed point de- 
scribes a smooth interface. But, since this fixed point 
is always unstable, the interface is always rough below 
d c = 2. (3) At the critical dimension the strong-coupling 
fixed point, as obtained from the two-loop calculation, 
tends to infinity. It is not clear, whether this divergence 
of the strong-coupling fixed point is just an artefact of 
the two-loop calculation or a general feature of any fi- 
nite order in perturbation theory. We suppose that it 
indicates that there is non-perturbative strong-coupling 
behavior for all d> 2. 

We have also highlighted the importance of carefully 
considering the structure of the dimensional regulariza- 
tion method. As explained in detail in Section III, some 
care is required in using this method, especially in dis- 
tinguishing between the dimensional dependence con- 
tained in 1/e poles and d-dependent factors originat- 
ing from purely geometric properties of the loop inte- 
grals. Whereas the 1/e poles constitute essentially noth- 
ing but a quite convenient way of keeping track of the 
UV divergences in the perturbation expansion, all other 
d-dependent factors characterize the symmetry and in- 
ternal structure of the model. Note that with these pre- 
cautions the results for the renormalization factors quite 
naturally obey the FDT in d = 1 dimension; it is not 
necessary to perform a separate treatment of this case. 

The two-loop calculations allow for no decisive conclu- 
sion about the existence of an upper critical dimension, 
at which mean field values for the exponents are recov- 
ered. To decide upon the existence of an upper critical 
dimension it would be highly useful to do a systematic 
expansion in 1/d within the framework of the Burgers- 
KPZ equation. To our knowledge, such an expansion is 
not yet available at the present time. 

Very recent results by Tu f38|, who has solved the 
mode-coupling equations without assuming any specific 
line shape, indicate that z < 2 for all dimensions d. 
However, the results obtained from mode-coupling the- 
ories are not quite conclusive, since it is not known in 
what sense such theories constitute a controlled expan- 
sion. Nevertheless, resummation techniques such as the 
mode-coupling theory are at present the only available 
analytical technique to study strong-coupling behavior. 

For the (1 + 1) -dimensional case there is very close 
agreement between numerical simulations and the results 
obtained from mode-coupling theory [p4| . Our two-loop 
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calculations demonstrate that for d = 1 all two-loop cor- 
rections vanish. This is a highly non-trivial result, be- 
cause for this to be the case, the vertex corrections enter- 
ing the calculations of the correlation and response func- 
tions must cancel each other. It is exactly those vertex 
corrections, which are neglected by mode-coupling the- 
ory, however. Hence our RG analysis provides a distinct 
hint why this approach has been so successful in 1 + 1 
dimension. Yet, further analysis is required to clarify this 
point @. 

Summarizing, it would thus be of interest to further in- 
vestigate the behavior of the strong-coupling fixed point 
by a systematic 1/d expansion, as well as by mode- 
coupling theory in high dimensions. In addition, it would 
be interesting to know how the value of the critical fixed 
point (obtained within the two-loop calculation) — es- 
pecially its d dependence — and the crossover scaling 
analysis in the smooth phase compare with results from 
numerical simulations. 
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V. APPENDIX 



In the appendices, we list the explicit results for the perturbation theory to two-loop order. Appendix A comprises 
the Feynman diagrams and the corresponding analytical expressions, while we present the evaluation of the relevant 
integrals in dimensional regularization in Appendix B. 



A. Two— loop Perturbation Theory for the Two— point Vertex Functions 

This appendix comprises the Feynman diagrams to two-loop order for the Burgers - Kardar-Parisi-Zhang equa- 
tion, and the corresponding momentum integrals. The integrations over the internal frequencies have already been 
performed using the residue theorem. 

We start with a list of the contributions to two-loop order to the fully wavevector- and frequency-dependent 
two-point vertex functions. In writing down the diagrammatic expansion for the dynamic functional one has to 
take into account restrictions which follow from causality. In section II we have not written down the Jacobian 
J[h] = T>[rj\/'D[h]. In general the Jacobian depends on the discretization of the Langevin equation (needed to give 
a proper definition to the path integral). As can be shown quite generally P4I the Jacobian cancels the equal time 
contractions of the field h and the response field h. Keeping in mind (or by choosing a discretization with the Jacobian 
equal to 1) one can omit the Jacobian in the dynamic functional. The Feynman diagrams, which account for the 
restrictions imposed by causality pSfl , are depicted in Fig. 7 and Fig. 8 for r^(q, cj) and T^ fi (q, w), respectively. 
In case of the former, adding the terms corresponding to both opposite time directions, respectively, provides some 
simplifications. The corresponding analytical expressions are: 



(a) + (b) 
(c) + (d) -- 



-2£> - ^iRe 



2ul 



(<z74-p 2 ) 



2\2 



\ 4 nD? 



-Re 



p (q/2 + p) 2 (q/2 - p) 2 iu/2v + q 2 /i + p 2 ' 
(q 2 /4-p 2 ) 2 f (q/4- p/2) 2 -fc 2 



8v$ J p (q/2 + p) 2 (q/2 - p) 4 J k iuj/2v + q 2 /4 + p 



(q/2-p) 2 [(q/4-p/2) 



2[(qk)/2-(pk)] 2 



(q/4 - p/2 + k) 2 (q/4 - p/2 - k) 2 [3(q/2 - p) 2 /4 + fc 2 ] 



2 + 



(q/2 - p)' 



iu/2u Q + (q/2 + p) 2 /2 + (q/2 - p) 2 /4 + fc 2 



1 



3(q/2-p) 2 /4 + fc ; 
iuj/2is + q 2 /4 + p 2 



(e) 



Re 



{q 2 /A-p 2 f 



1 



(q/2 + p) 2 (q/2 - p) 2 ilu/2is q + q 2 /A + p 2 

[(q/4- p/2) 2 -fc 2 ] 2 

k (q/4 - p/2 + k) 2 (q/4 - p/2 - k) 2 [3(q/2 - p) 2 /4 + fc 2 ] 



x 2 



(q/2 - p) 5 



(/KG?) — 



x Re 



A^o 3 



Re 
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q 2 /A - (p + k) 2 
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1 1 (q/2 - p) 2 



(h) + (<) 



x Re 



'-Re 



q 2 /4-p 2 



p (q/2 + p) 2 (q/2-p) 



(qk) 2 /4-(pk) 2 

fc 2 



g 2 /4-(p + k) 



1 
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iuj/2v + (q/2 + p) 2 /2 + (q/2 - p) 2 /4 + k 2 J ' 
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(5.1) 
(5.2) 



(5.3) 



(5.4) 



icu/2is + q 2 /A + p 2 iLu/2v + q 2 /A + p 2 + (q/2 - p - k) 2 + (q/2 + p) 2 + k 2 



(5.5) 



iu/2v + q 2 /4 + p 2 iuj/2v + q 2 /4 + p 2 iu/v + (q/2 - p - k) 2 + (q/2 + p) 2 + fc 2 
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(c) 



A^o 2 A 9 2 /4-p 2 9 2 + 2(qp) /■ (q/4 -p/2) 2 



32^ 5 7 p (q/2 + p) 2 (iu/2v + 9 2 /4 + P 2 ) 2 h (q/4 - p/2 + k) 2 (q/4 - p/2 - k) 2 
(q/2 - p) 2 [(q/4 - p/2) 2 + fc 2 ] - 2[(qk)/2 - (pk)] 2 
iu;/2v + (q/2 + p) 2 /2 + (q/2 - p) 2 /4 + fc 2 ' 



AqDq /■ g 2 /4-p 2 g 2 -2(qp) 
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iiv/2v + <7 2 /4 + (p + k) 2 iu/vo + (q/2 - p - k) 2 + (q/2 + p) 2 + fc 2 



m _ KDp f g 2 + 2(q P ) f (qk) 2 /4-(pk) 2 

U 8^ 5 J p (q/2 + p) 2 (iw/2n, + 9 2 /4 + p 2 ) h fc 2 

g 2 /4-(p + k) 2 1 



ioj/2v + 9 2 /4 + (p + k) 2 iu/uQ + (q/2 - p - k) 2 + (q/2 + p) 2 + fc 2 
q 2 - 2(qp) /■ [(qk) 2 /4 - (pk) 2 ] [g 2 /4 - (p + k) 2 ] 



Aq-Pq y g 2 - 2( qP ) r 

8^ 5 j p (q/2 - p) 2 ^/2^ + 9 2 /4 + p 2 ) h 



fc 2 [iw/2^ + ( j 2 /4+(p + k) 2 ] 



( 1+ 2(q/2-p) 2 



(q/2 - p - k) 2 + (q/2 - p) 2 + fc 2 V iu/v + (q/2 - p - k) 2 + (q/2 + p) 2 + fc 2 , 
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(k) = 



x 4 n 2 



9 2 -2(qp) f (qk)/2 + (pk) 



g 2 /4-(p + k) 2 



4i/« 7 p (q/2 - p) 2 j fc fc 2 (q/2 p k) 2 ^/2^ + <Z 2 /4 + (p + k) 2 

(gk)/2 - (pk) - fc 2 (q/2 - p) 2 

(q/2 - p - k) 2 + (q/2 - p) 2 + fc 2 iu/v Q + (q/2 - p - k) 2 + (q/2 + p) 2 + fc 2 ' 



(5.15) 



Upon collecting the several contributions, these expressions simplify considerably. One finds that the vertex 
functions may be split into UV-singular and UV-regular parts according to = r~| 9 + r~ 4 ™ 9 , and similarly 



hh 



^kh = ^}i h 9 + ^r™ 9 - After some tedious, but elementary algebra one arrives at the following explicit results: 



hh 



r^(q,^r 



Re 



(iu + v$q 2 ) 



q 2 /A-p 2 



Ik i 



1 



4^ 7 J p (q/2 + p) 2 ilu/2v + q 2 /A + p 2 
(qk)/2 + (pk) „ 



iu/2u + q 2 /A+(p + k) 2 fc 2 (q/2 - p - k) 2 

(qk)/2 - (pk) - fc 2 



+ 



q 2 /A - p 2 



iu/vo + (q/2 - p - k) 2 + (q/2 + p) 2 + k 2 
1 (qk)/2 + (pk) 



(q/2 + p) 2 (q/2 - p) 2 J k iu/2v + q 2 /A + (p + k) 2 fc 2 (q/2 p - k) 2 



Re- 



(q/2 - p) 2 fc 2 - [(qk)/2 - (pk)] s 

(q/2 - p - k) 2 + (q/2 - p) 2 + k 2 V^iu/2v + <? 2 /4 + p 2 
1 (q/2 - p) 



+ 



+ 



?hhM Sin9 
A^o 3 



-2D a 



ilu/2is + q 2 /A + p 2 iu/uo + (q/2 - p - k) 2 + (q/2 + p) 2 + k 2 
(q 2 /A-p 2 ) 2 1 



\ 2 n 2 
^Rc 
2v% 



A^l 



■Re 



(q 2 /A-p 2 ) 2 



p (q/2 + p) 2 (q/2 - p) 2 iu/2v + q 2 /A + p 2 

(q/4 - p/2) 2 - k 2 



Re 



(q/2 + p) 2 (iu/2v + q 2 /A + p 2 ) 2 J k (q/4 - p/2 + k) 2 (q/4 - p/2 - k) 2 
(<Z 2 /4-p 2 ) 2 1 



(q/2 + p) 2 (q/2 - p) 2 (iu/2u + <? 2 /4 + p 2 ) 2 

(q/2-p) 2 fc 2 -[(qk)/2-(pk)] 2 



■^Re 



(<Z 2 /4^P 2 ) 2 



k (q/4-p/2 + k) 2 (q/4-p/2-k) 2 
1 



p (q/2 + p) 2 (q/2 - p)* ilo/2is + q 2 /A + p 2 

(q/2 - p) 2 k 2 - [(qk)/2 - (pk)] 2 



r^(q,^) r 



-(iu + v q 2 ) 



(q/4-p/2 + k) 2 (q/4-p/2-k) 2 
g 2 (g 2 /4 + p 2 ) - 2(qp) 2 
' p (q/2 + p) 2 ^/2^o + g 2 /4 + p 2 ) 
(qk)/2 + (pk) „ 



k iuj/2v + rj 2 /4 + (p + k) 2 fc 2 (q/2 - p - k) 2 

(qk)/2 - (pk) - k 2 



' iu/vo + (q/2 - p - k) 2 + (q/2 + p) 2 + k 2 



+ 



x 4 r> 2 



g 2 (g 2 /4 + p 2 )-2(qp) 2 



Av% J p (q/2 + p) 2 (q/2 - p) 2 ^/2^ + <7 2 /4 + p 2 ) J k ilu/2v + <? 2 /4 + (p + k) 2 
(qk)/2 + (pk) (q/2 - p) 2 fc 2 - [(qk)/2 - (pk)] 2 

fc 2 (q/2 - p - k) 2 ilu/vo + (q/2 - p - k) 2 + (q/2 + p) 2 + k 2 
XiD 2 fq 2 -2(qp) f 1 (qk)/2 + (pk) ^ 



KDq [ q- 2(qp) f 
Av% J p (q/2 - p) 2 J k i 



iu/2v + q 2 /A + (p + k) 2 fc 2 (q/2 - p - k) 2 
(q/2 - p) 2 fc 2 - [(qk)/2 - (pk)] 2 



[(q/2 - p - k) 2 + (q/2 - p) 2 + fc 2 ] [iu/u + (q/2 - p - k) 2 + (q/2 + p) 2 + fc 2 ] 



(5.16) 



(5.17) 



(5.18) 
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\)^0 



q 2 /4 - p 2 



q 2 {q 2 /A + p 2 ) - 2(qp)^ 



Av 2 J p (q/2 + p) 2 (q/2-p) 2 iu/2v + q 2 /4 + p 2 



(5.19) 



XtD 2 f q 2 /4-p 2 q 2 (q 2 /4 + p 2 )-2(qp) 



16^ J p (q/2 + p) 2 (iw/2m + g 2 /4 + p 2 ) 2 J fe (q/4 - p/2 + k) 2 (q/4 - p/2 - k) 2 



(q/4-p/2) 2 



x 4 r> 2 



g 2 /4 - p 2 



q 2 (q 2 /4 + p 2 )-2(c l py 



8i/« J p (q/2 + p) 2 (q/2-p) 2 (tw/21* + g 2 /4 + p 2 ) 2 ' 

(q/2 - p) 2 fc 2 - [(qk)/2 - (pk)] 2 



A^ 2 f q 2 /4-p 2 



2(qp) 



(q/4-p/2 + k) 2 (q/4-p/2-k) 2 

(q/2 - p) 2 fc 2 - [(qk)/2 - (pk)] 2 



(q/2 - p)4 iu/2vo + q 2 /A+p 2 J k (q/4 - p/2 + k) 2 (q/4 - p/2 - k) 2 
q 2 /A _ p 2 q 2_ 2(qp) 



16^ 5 Jp (q/2 " P) 2 iw/2wb + g 2 /4 + P 2 



(q/2 - p) 2 fc 2 - [(qk)/2 - (pk)] 5 



k (q/4 - p/2 + k) 2 (q/4 - p/2 - k) 2 [3(q/2 - p) 2 /4 + k 2 ] 



A powerful check to these lengthy calculations is the investigation of the situation at d = 1, where the above 
expressions simplify considerably. Using some elementary algebra again, one finds at d = 1 



voq 



(5.20) 



which ensures the validity of the fluctuation-dissipation theorem, for both renormalized mom entum- and frequency- 
dependent quantities v(q,uj) and D(q,u>), defined as the left- and right-hand sides of Eq. ( 5.20] ), respectively, coincide. 
The explicit result reads: 



v(q,u)= is 



A 2 , A) 

41/g 



Re 



An-Dn 



p ^/2^ + g 2 /4 + p 2 

(g/2-p) 2 



(5.21) 



16^ 6 



Re 



(iw/2j/ + g 2 /4 + p 2 ) 2 J k k(q/2 -p-k) 



\ 4 n 2 
8i/, 



JiRejW^ 2 ) /t 
o V Jp ! 



g/2-p 



itv/2v + q 2 /4 + p 2 



1 



fc zcj/2i/ + <z 2 /4 + (p + A;) 2 
1 



q/2- p-k iu/va + (q/2 - p - k) 2 + ( ? /2 + p) 2 + fc 2 

which already is a useful and interesting result of the two-loop calculation on its own standing. 

We now return to the general d-dimensional case. In evaluating the UV-singular contributions, one has to be 
careful to choose a normalization point (NP) where either q or to are finite, in order not to interfere with the IR 
singularities, which would show up as poles in l/(d — 2), too. A convenient choice is NP: q = 0, iuj/2v = /i 2 ; thus we 
find 



r^(q^) 



sing 



NP 



-2D 



l+3o 



+23o 2 



+43o 



1 



p fi 2 Z„+p 2 
1 



+ 3o 2 



P 



p 2 /A - k 2 



P (p 2 +P 2 ?J k (p/2 + k) 2 (p/2-k) 2 
p 2 k 2 - (pk) 2 



(V 2 +P 2 ) 2 J k (p/2 + k) 2 (p/2-k) 2 



p 2 k 2 - (pk)' 



p p 2 (fi 2 +p 2 ) 



(p/2 + k) 2 (p/2-k) 5 



0(gl) 



(5.22) 



Ar 

dq 2 hh 



sing 



NP 



d-2 



d-2 



-go 



p n 2 Z v +p 2 
p 2 /A - k 2 



p(v 2 +P 2 ) 2 J k (p/2 + k) 2 (p/2-k) 5 
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d-2 2 



1 



p 2 k 2 - (pk) z 



d-2 n 
-2- - 2 



P {p 2 +P 2 ) 2 Jk (p/2 + k) 2 (p/2-k) 2 

i r P 2 k 2 - ( P k) 2 



-9l 



P P 2 ^ 2 +P 2 )Jk (p/2 + k) 2 (p/2-k) 2 

i r P 2 k 2 - ( P k) 2 



P M 2 + P 2 J k (p/2 + k) 2 (p/2 - k) 2 (3p 2 /4 + k 2 ) 



* 2 
-d 9 ° 

4 2 



1 



[p 2 fc 2 - (pk) 2 ]fc 2 



p V 2 +P 2 Jk (P/2 + k) 2 (p/2 - k) 2 (3p 2 /4 + fc 2 ) 2 
1 /" (p 2 /4- fc 2 )[p 2 fc 2 - (pk) 2 ](pk) 2 



p2( M 2 +p 2) J k (p/2 + k) 4 (p/2 - k) 4 (3p 2 /4 + k 2 ) 



O(gl) 



For the derivation of the latter expression, the following relations have proven very useful: 



(qp)/(p)= 



(qp)(qk)(pk)/(p,k) = 



P 2 .f(p) ■ 



p J k 



p J k 



( P k) 2 /(p,k). 



(5.23) 

(5.24) 
(5.25) 



B. Integrals in dimensional regularization 



In this appendix, we list the results for the above integrals, as obtained from the dimensional regularization scheme 
23151. 



1 



P M 2 +P 2 
P 2 



p 2 /A - k 2 



p (fi 2 +P 2 ) 2 J k (p/2 + k) 2 (p/2-k) 2 



id -I) 



C 2 ^T(l-e)T(l + e/2) 



1 



p 2 k 2 - (pk) 2 



(H 2 +P 2 ) 2 J k (p/2 + k) 2 (p/2-k)- 



= (d-l) 



1 



p 2 k 2 - (pk) 5 



p 2 {(J. 2 +P 2 ) J k (p/2 + k) 2 (p/2-k)- 
p 2 fc 2 - (pk) 2 



= (d-l) 



2e 2 r(l-e/2) 

c 2 ^r(i- £ )r(i + e /2) 

4e 2 r(l-e/2) 
C 2 ^ r(l-e)r(l + e/2) 



4e 2 (l + e) r(l-e/2) 



M 2 +P 2 Jk (p/2 + k) 2 (p/2-k) 2 (3p 2 /4 + fc 2 ) 



c 2 ^ r(i- c )r(i + £ ) 
8 £ r(i- e /2)r(i + e /2) /oo(d) ' 



i 



[p 2 fc 2 - (pk) 2 ]fc 2 



(p/2 + k) 2 (p/2 - k) 2 (3p 2 /4 + fc 2 ) 2 

c 2 ^ r(i-e)r(i + c ) 



M 2 + P 2 Jk 
= -(d-l) 
1 

p 2 (M 2 +p 2 ) 7 fe (p/2 + k) 4 (p/2 - k) 4 (3p 2 /4 + fc 2 ) 
(2-e)C 2 ^ r(l-e)r(l + e) 



16e r(l - e/2)r(l + e/2) 

(p 2 /4 - k 2 )[p 2 k 2 - (pk) 2 ](pk) 2 



(d + e)I w (d) + ^hi(d) + ^ /21(d) 



(d-1)- 



32e r(l - e/2)r(l + e/2) 



(3-e)/oo(d)-i5_J!/ u ( d ) 



3 e / j\ 7(4 - e) _ 

/21(d) + -J. 

4 u ' 16 



22(d) -i^ijg^ + i 



x(l — x) 



■^[i 01 (d)-i 12 (d)} 



+ (d-i) 



(l + rr-a; 2 ) 2 ^/ 2 

c d V £ r(i- e )r(i + e /2) / 2 



d.r 



16e 2 



r(l-e/2) 



1 + e 



(2-e) 



(5.26) 
(5.27) 
(5.28) 
(5.29) 

(5.30) 
(5.31) 



(5.32) 
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Here C d = T(2 - d/2)/2 d - 1 w d / 2 is a geometry factor, while 

' '■' v(i-y) r 



M) = I I 

Jo Jo 



[x(l - x)y 2 + (1 - y)(3 + y )/4\s+2-d/2 
and 



dxdy 



I rs (d) 



IrM ~l l [ X (l-x)+7 r yY + ^ dxdy 



are parameter integrals emerging upon the use of Feynman parametrization. Finally, we note that 

r(l-e)r(l + e/2) = 2 
r(l-e/2) L + U ^ e h 

r(i- e )r(i + e) 



r(l-c/2)r(l + e/2) 

and 



/ (l+x-x 2 ) 2 _ 5^5 y/l-l 5 

are to be used when performing the e expansion leading to the results for the Z factors. Note also that at d 

87 00 (2) - 47 10 (2) - 217 n (2) + J 21 (2) + 77 22 (2) - 27 32 (2) 

-47 01 (2) + 47 12 (2) + In ^i±± -1+4 = 8- F„(2) = . 
5V5 v5 - 1 5 
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Figure Captions: 



FIG.l: Basic elements of the dynamical perturbation theory for the Burgers - Kardar-Parisi-Zhang equation, (a) 
Correlation and response propagators, (b) Three-point vertex. 



FIG. 2: One-loop fixed point of the Burgers-KPZ equation as function of the surface dimension d. The divergence 
at d — 3/2 turns out to be an artefact of the one-loop approximation. 



FIG. 3: Two-loop fixed points of the Burgers-KPZ equation as functions of d (full lines). It is not clear if the 
divergences at d = 1 and d = 2 are inherent to the low order of perturbation expansion used, as was the case for the 
divergence of the one-loop fixed point (dashed here) at d = 3/2, or if they rather represent a generic feature of the 
model. 



FIG. 4: The correlation length £ in the smooth phase versus \g — g*\/g'c- The correlation length crosses over from 
£ oc g 1 ^ at small couplings to £ oc \g — g*\~ v as g approaches the critical coupling g*. In the figure we have set e = 1. 



FIG. 5: The renormalized noise amplitude D{g) (solid line) and surface tension u(g) (dashed line) (up to a nonuni- 
versal amplitude) in the smooth phase versus \g — g* \/g*. The noise amplitude shows a crossover from an exponential 
increase D(g) oc e 9 ^ e at small values of g to a power law divergence D(g) oc \g — g* as g approaches the critical 
coupling g*. The renormalized surface tension v(g) crosses over from exponentially decreasing v(g) oc e~ g l 2 to a 
constant at criticality. The curves are plotted for d = 3 (i.e. e = 1). 



FIG. 6: The characteristic frequency u> c (g) ( U P to a nonuniversal amplitude) in the smooth phase versus \g — g*\/g*- 
The characteristic frequency crosses over from u c (g) oc g~ 2 / e e~ 9 / 2 at small couplings to u c (g) oc \g — g*\~ 2 as g 
approaches the critical coupling g*. The figure shows the result for e = 1. 



FIG. 7: Feynman diagrams for the dynamical perturbation expansion of r^(q, u>) to two-loop order 



FIG. 8: Feynman diagrams for the dynamical perturbation expansion of r^(q, u>) to two-loop order 
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